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Abstract 

This paper addresses the problem of designing suitable 
satellite constellations for Radio Occupation missions. Radio 
occultation for the Earth atmosphere usually requires global 
coverage and short interval measurements. The 2-D Lattice 
Flower Constellations theory is here applied to design 
constellations in order to maximize active time by providing 
global coverage with frequent measurements. Optimizations 
are performed using Genetic Algorithms to estimate 
constellation design parameters. Optimization is constrained 
by altitude range (to limit atmospheric drag and avoid Van 
Allen belts) while orbit nodal precession is used to obtain 
global coverage. The resulting constellation geometries have 
been explored from coverage performance perspective. Two 
coverage estimation methods have been used to evaluate 
coverage performances. Performances of some solutions are 
provided to make optimality selection. 
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Introduction 

Radio Occultation (RO) is a recently proposed remote 
sensing technique to measure, in real time, physical 
parameters of the atmosphere, such as density and 
water vapor. The technique requires at least two 
satellites, one transmitting an electromagnetic signal 
and the other receiving it. The signal arriving to the 
receiver is refracted as it travels through the 
atmosphere, and the refraction magnitude is a 
function of the density and water vapor. Detailed 
information on radio occultation can be found in Ref. 
[Yunck 2002]. Since the GPS satellites have already 
provided this electromagnetic (radio) signal, 
techniques known as GPSRO are proposed for LEO 
satellites [Kursinski et al. 1997,Melbourneet al. 1994]. 
For GPSRO missions, the signal bending (due to the 
atmosphere refraction) is measured by detecting 
Doppler shift. In particular, it is possible to retrieve 
information on temperature and pressure using an 



Abel transform [Bracewell 1965] for atmosphere 
altitudes below the ionosphere. 

The name, radio occultation, comes from the odds for 
two satellites to see each other just before they are 
occulted by the Earth during the relative orbital motion. 
In that short time, when they can see each other 
through the atmosphere, the electromagnetic signal is 
sent and the radio occultation measurement can be 
performed. In general, for two satellites having 
different orbital periods, the RO measurement time 
happens only once per orbit. In this paper, we propose 
three different satellite constellations where a subset of 
satellite pairs sees each other for almost all the orbital 
period. In particular, these constellations, using 
circular orbits only, are analyzed and compared in 
terms of atmosphere coverage without consideration 
on the atmospheric refraction. When the actual 
atmospheric refraction is taken into account, the 
radius of the optimal orbit is slightly changed but the 
atmospheric coverage remains substantially identical. 

When a radio signal from a transmitter passes through 
the atmosphere, its phase is perturbed in a way related 
to the atmospheric refractivity along the 
electromagnetic signal path. In general, the refraction 
angle is obtained from phase measurements on a 
receiver. Then, atmospheric quantities such as density, 
temperature, pressure, and moisture can be derived 
from it. 

Radio Occultation missions can be classified into the 
following three categories: 

1. two (or more) satellites on the same orbit (forming 
"string-of-pearls"), 

2. geostationary satellites with Low Earth Orbit 
(LEO)satellites, and 

3. Global Navigation Satellite System (GNSS) with 
LEOsatellites [Liou et al. 2010]. 

The first RO experiment was performed during the 
occultation of the Mariner V spacecraft by Venus on 19 
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October 1967 after theoretical studies [Fjeldbo et al. 
1965, 1971]. RO investigation using geostationary 
satellite with LEO satellite was made by radio link 
between two satellites in 1974 [Liu 1978]. The first was 
the Applications Technology Satellite (ATS-6) and the 
second was Geodetic Earth Orbiting Satellite (GEOS-3). 
Because these approaches provide limited data, many 
satellites are required for Earth's atmosphere global 
coverage, leading to unaffordable costs. This situation 
changed when Global Positioning System (GPS) 
signals were proposed to be used for RO. The RO 
technique has been applied to GPS signals as observed 
by LEO satellites in order to monitor the Earth's 
atmosphere with acceptable costs. The first experiment 
of GNSS-LEO was the GPS/MET (METeorology) 
experiment using the Microlab-1 satellite in 1995 
[Melbourne et al. 1994, Ware et al. 1996]. Then, several 
other RO experiments have been conducted using GPS 
signals, notably CHAMP (Challenging Minisatellite 
Payload), GRACE (Gravity Recovery And Climate 
Experiment), and the FORMOSAT-3/COSMIC 
(Constellation Observing System for Meteorology, 
Ionosphere, and Climate) constellation [Wickert et al. 
2001, Wickert et al. 2005]. Then, several RO missions 
using GNSS-LEO technique have been designed for 
the Russian GLONASS and the European Galileo 
constellations. With the advent of more efficient 
technologies, the use of cost-effective small satellites 
allows LEO-to-LEO links to provide global coverage 
with acceptable mission costs. 

Optimal constellation design is a difficult problem 
because of the high dimensionality of the problem as 
each orbit depends on six parameters. However, it is 
possible to design different configurations of LEO 
satellites (LEO constellations) able to obtain high RO 
active times and global coverage. The optimal design 
of such constellations can be done using existing, well- 
proven, satellite constellation frameworks (e.g., 
Streets-of-coverage [Liiders 1961, Liiders et al. 1974] 
and Walker's [Walker 1971, 1974]) as well as the 
recently proposed 2-D Lattice Flower Constellations 
(LFC) [Avendano et al. 2010]. The goal is to obtain a 
reasonable reduced problem dimensionality without 
cutting out potential useful solutions. For elliptical 
constellations, the LFC requires seven parameters 
three of which are integers; while for circular 
constellations, there are there integers in five 
parameters because eccentricity and argument of 
perigee are not necessary. The optimization of LFC 
implies that these five parameters have to be found to 
fulfill mission requirements. This research study 



shows how to derive these parameters in order to 
maximize the atmosphere coverage. Reasonable initial 
radius ranges can be obtained by limiting the 
atmospheric drag and avoiding the inner Van Allen 
belt. 

This paper is organized as follows. The first section of 
this paper briefly describes the radio occultation 
mission. Then, a summary on the 2-D LFC theory is 
provided. Then design considerations and parameters 
for LFC and optimize the design of the constellations 
are discussed. Finally, theperformance of three 
selected constellations for radio occultation mission 
are presented. 

Radio Occultation Mission 

Mission Needs 

The cost function defining the optimal RO 
constellation for global coverage of the Earth's 
atmosphere requires the maximization of the odds at 
which pairs of satellites can see each other through the 
atmosphere (active or observation time). In the next 
section, the basic geometry of two satellites is defined 
during the RO active time. 

Geometry of Radio Occultation 




FIG. 1 TWO SATELLITE GEOMETRY DURING RADIO 
OCCULTATION ACTIVE TIME 

The simplest optimal RO constellation is made of 
satellites on the same circular orbit of radius R. Fig.l 
shows how the maximum central angle (ctmax) between 
two satellites for RO mission depends on the 
minimum atmosphere radius (Rmin) while minimum 
central angle (amm) depends on maximum atmosphere 
radius (Rmax). In this study we consider the atmosphere 
altitude ranging from 20 km to 200 km. Two satellites 
on the same circular orbit see each other through 
atmosphere if 
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Then, using Eq. (3), we obtain 



where a „ = 2 cos" 1 1 ] , a_ = 2 cos" 1 f I and 



R 



R 



R. and R . are the position vectors of the i -th and the 
j -th satellite, respectively. 

Atmosphere Coverage 

1) Linear Swath Geometry of Two Satellites at a 
Specific Time 

Linear swath between two satellites at a given 
instant time is as shown in Fig. 2. The satellites are 
located at orbital positions S, and S . , while the 

Earth's center is indicated by the point O . The 
distance between the two satellites is L and the 
angles a and a* are the central angles ZS,OS / and 

ZS*OS* , respectively. 




FIG. 2 LINEAR SWATH GEOMETRY FOR TWO SATELLITES AT A 
GIVEN TIME 

To analyze the coverage, let us assume the radio 
signalis a line segment passing through the 
atmosphere with norefraction (which is practically 
very small). The end pointsof this segment are 
indicated by S, and S y , respectively .Midpoint of 
the segment is point M . The length / of the linear 
swath is simply given by 



: 2R... 



1 ' 

1-COS 











(2) 



where Rmax is maximum atmosphere height. 
Considering the two triangles, AOMSj and 
AOMS*, we can write AOMS*, we can write 



OM = 7? cos 



2 } = *-«» 



(3) 



l = 2RM- 



R (a 

5 



■cos 



2.R1.-R 



cos a + 1 



(4) 



where R is orbital radius. 



2) Planar Swath Geometry of Two Satellites 
During a Time Interval 

Fig. 3 shows planar swath geometry for two 
satellitesduring a time step (t k+1 -t k ) . Considering 
two satellitesin time t k and t k+1 , we have 4 position 

vectors of satellites [S t (t k ), Sj(t k ), Sj(t k+1 ), S ; (i w )]. 




FIG. 3 PLANAR SWATH GEOMETRY FOR TWO SATELLITES 
DURING A TIME STEP 

These position vectors indicate 4 vertices of planar 
swath that can be obtained as follows. 



s;(* t )=s ( (t t ) 

S%) = S j (t k ) 
S](t k+1 ) = S j (t k+1 )- 

s;(t t+ i)=s l (» l+I )- 



Sj(t k )-S,(t t ) L-l 

"is.a t )-s,.(^)i 2 

S,.(f t )-S.(f t ) L -i 
IS.(g-S.(? t )l 2 

s,.(f t+1 )-s.fc +1 ) L'-r 

IS^-S.C^)! 2 



(5) 



IS y (? t+1 )-S,.(f, +1 )l 2 

These 4 points, [S*(t k ), S)(t k ), S)(t M ), S](t k+l )] , 

identify a quadrilateral planar swath. The time step 
is small enough to assume that satellites move in 
straight line. Since this quadrilateral is skew, these 
points are transformed to topocentric-horizontal 
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coordinate system with projection to surface. 

3) First Coverage Estimation-sum of Areas of 
Simple and Complex Quadrilaterals 

There are two quadrilateral cases as shown in Fig. 4. 
The first one is simple quadrilateral while the 
other is a complex (self -intersecting) quadrilateral. 



s> t ) s 



v, 

down & 



h, = U,xU 2 
— ®down 



h, = V.xV, 




V. 

FIG. 4 SIMPLE AND COMPLEX QUADRILATERALS 

Using position vectors which are given by Eq. (5) 
we can get vectors such as Ui, U2, Vi and V2. Then, 
vectors hi and I12 can be obtained by cross products 
as shown in Fig. 4. The condition associated with a 
simple quadrilateral is 



h 2 >0 



(6) 



Because all quadrilaterals consist of two triangles 
regardless of shape, area of simple or complex 
quadrilaterals can always be obtained by summing 
the areas of two triangles. Area for simple case is 
obtained by 



J h, I I h J 

Simple - ~ 



(7) 



In complex case, there are some additional 
processes because we find intersecting point (Pi) 
and then define related vectors as follows. 



u;=p,- s*(rj 
v- Sim p, 
h; = u, x u' 2 

Similarly, area for complex case is given by 
I h! I I h' I 

A — 1 j r_ 

complex ~ ^ 



(8) 



(9) 



The total coverage is obtained by summing all 
planar swaths by all satellites pairs combinations. 
For instance, by summing all planar swaths by all 



satellite combinations during one quart of the 
period gives a quarter period coverage. However, 
net coverage cannot be obtained this way because 
of the existing overlapping areas. In order to 
estimate the net coverage, methods to count 
overlapping area are required. 

4) Second Coverage Estimation - Using Uniform 
Dis tributed Points on a Sphere 

The creation of uniform distribution of points on a 
sphere can be used to estimate net coverage by 
checking the number of overlapping times. Among 
all the existing methods to generate uniformly 
distributed points on a sphere the method 
presented in Ref. [Mortari et al. 2011] is adopted 
here. Starting with an icosahedron (20 identical 
equilateral triangular faces) and performing 8 
sequential divisions in identical spherical triangles, 
20 -2 s =5,120 triangles with same spherical areas 
are obtained whose centers approximate 5,120 
uniformly distributed points on a sphere, as shown 
in Fig. 5. 




FIG. 5 ICOSAHEDRON SPLIT BY 8 DIVISIONS 

Because every planar swath is quadrilateral, it can 
be divided into two triangles. Therefore, checking 
if points are inside or outside triangle can 
determine if points are observed or not. Fig. 6 
shows the geometry of observation and non- 
observation cases. In case of observation (left), the 
sum of the areas of the three sub-triangles is equal 
to the area of the original triangle while the non- 
observation case (right) is experienced if the total 
area is greater than the area of the original triangle. 





s,v,) 



FIG. 6 GEOMETRY OF OBSERVATION AND NON- 
OBSERVATION CASES 



si, .;, 
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In order to check observation condition, all 
uniformly distributed points are transformed to 
topocentric-horizon coordinate system with mid- 
point of triangle as origin, and then projected to 
surface as all points of quadrilateral are 
transformed. Equation to check whether the point 
Fv is observed is given by. 



a(p v s'(t k ) s*a t ))+A(p v s*(fjs*(f t+1 )) 

+A(P v S*(? t+1 ) S*(f,)) = A(S*a t ) S](t k ) S](t k+1 )) 



(10) 



2-D Lattice Flower Constellations Design 

The 2-D Lattice Flower Constellations (LFC) design 
methodology has been introduced in Ref. [16]. In 
general LFC are characterized by four continuous 
parameters (semi-major axis, eccentricity, inclination 
and argument of perigee) and three independent 
integer parameters establishing the constellation 
satellite distribution in the (M,D.) -space. These integer 
parameters are: the number of orbital planes, N , the 
number of satellites per orbit, N so , and the 
configuration number, N c (phasing parameter). Using 

these integer parameters the satellites' right ascension 
of the ascending node (fi r ) and the initial mean 

anomaly (M 1} ) are solutions of the following equation 



No 
N, 







= In 



i-1 
7-1 



(11) 



where i = l,-,N B , j = l,-,N m , and N c e [l,N a ] . The 
"i-j" element is the j-th satellite on the i -th orbital 
plane. If repeating ground tracks are required, then 
the compatibility equation 

2x ™ 2x 



n + co 



- = N d T d =N d 



(12) 



where T p is the orbit nodal period and T d is the 

Greenwich nodal period, provides the value of the 
orbit radius for each coprime integers, N and N d . Eq. 

(12) takes into account the main gravitational 
perturbation due to the Earth oblateness, known as the 
J 2 effect. The secular and persistent J 2 effect modifies 
the mean motion according to 



n = n n 



\ + -J 7 

4 2 



(2-3sin 2 ^•)^/^V 



(13) 



where n = J-^- is the unperturbed mean motion, 
v a 

linearly changes the right ascension of the ascending 



node, 



n = — J, 



n cos i , 



as well as the argument of perigee, 



3 { R 

03 = — J 2 — | n(5cos 2 i-l). 

4 { p 



(14) 



(15) 



The LFC theory allows us to design constellations 
using circular or small eccentric orbits, within the 
range of validity of Carter's definition of compatibility, 
given in Eq. (12). For non small eccentric orbits, the J 2 
perturbation slowly destroys the constellation if no 
control is used to compensate the rotation of the apse 
line. However, in this paper circular orbits will be 
adopted, therefore, we can set e = while co (not 
defined for circular orbits) can be set to any constant 
value. Then, the number of continuous parameters for 
LFC is reduced from 4 to 2. In addition, all satellites in 
the constellation have same semi-major axis and 
inclination in order to avoid different orbital nodal 
regression. 

Minimum Distance Constraint 

To avoid the design of constellations with satellites 
colliding the results provided by Ref. [Speckman et al. 
1990] is adopted. In that analytical study, the closest 
approach between the two satellites (p min ) in two 
circular orbits with the same radius and inclination, is 
analytically expressed by the following equation 



1 + cos 2 i + sin 2 /cosAQ . f AF 

sin 

2 { 2 



(16) 



Where 

AF = AM -2 tan -1 



. ( AD. 

-cos (tan 

I 2 



and where AM and An are the differences in mean 
anomaly and right ascension of ascending node, 
respectively. Note that p min must be scaled by the orbit 
radius to find the actual value of the minimum 
approach distance. Due to the regular pattern (lattice) 
of the LFC, it is not necessary to evaluate the 
minimum distance using all pairs of satellites. It is 
sufficient to evaluate the minimum distance between 
the first satellite [Cl n ,M i J with all the other satellites 
staying on different orbital planes. This greatly 
simplifies the effort in the optimization to avoid 
constellations affected by satellite conjunctions 
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(sometime common for symmetric distribution). 
Constraint that no two satellites are ever closer than 
half the distance between two consecutive satellites in 
the same orbit has been used. For the specific case of 
N*, = 1 (just one satellite per orbit), the constraint 
adopted is that the minimum distance between any 
pair of satellites is greater than 100 km. 

LFC Optimization 

1) Genetic Algorithms 

Genetic Algorithms are empirical methods to 
solving optimization problems by mimic the 
natural selection/mutation process driving the 
biological evolution. More information about 
Genetic Algorithms is available in Ref. [Goldberg 
1989]. Although there is no guarantee that Genetic 
Algorithms provide the optimal solution (and this 
is true for most optimization methods), Genetic 
Algorithms are most suited in multi-parameter and 
highly nonlinear problems [Gallagher et al. 1994]. 
Because constellation design for radio occultation is 
a highly nonlinear problem, a Genetic Algorithm 
has been selected to find the optimal constellation 
parameters for the given set of mission 
requirements. 

2) Fitness Function 

The fitness function utilized to drive the 
optimization process is designed to maximize the 
percentage of active time with respect to a 
reference time. Fitness function has been 
implemented to count the times that pairs of 
satellites see each other through the atmosphere. 
This is simply determined by the geometric 
conditions expressed by Eq. (1). For computation of 
the percentage of active time, geometric conditions 
for all the i-j satellite combinations should be 
considered. Thus, the optimality is defined by the 
minimization of the following function 

L.ffi (17, 

where K is the number of times satisfying 
geometric conditions, Afis the time step, and ris 
the orbital period. 

3) Design Parameter Constraints 

In order to run the optimization process the design 
space should be first defined. Acceptable orbit 
radius is bounded by the maximum atmosphere 
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altitude (to limit the atmospheric drag) and by the 
inner Van Allen radiation belt. The maximum 
altitude is set to h msa = 1, 000 km because the lower 
boundary of the inner Van Allen Radiation belt is 
located at an altitude of about 1,000 km [Gilmore 
2002]. Considering atmospheric drag, the 
minimum altitude is set to h^ n = 300 km. Lower 
altitude are impractical due to the larger 
atmospheric drag [Young et al. 2000]. Using these 
two altitude constraints the maximum and 
minimum periods are 

Tmm = 2;r ^-y and T _ = 2}T yy , 

where and fi are the Earth's radius and 
gravitational constant, respectively. The ranges of 
these parameters are given in Table 1. 



TABLE 1 PARAMETERS RANGE 



Parameters 


Range 


[N„,N S0 ] 


[1,12] [2,6] [3,4] [4,3] [6,2] [12,1] 


N c 


[1, N a ] 


T 


[1.51 hr, 1.75 hr] 


i 


[0, 180 deg] 



Re-oriented Constellation 

Any constellation can be re-oriented just by using 
rigid rotation matrix with free selection of axis and 
angle of rotation in order to provide desired 
performance in specific region. In Fig. 7, n^, , e k and h t 
are node, eccentricity, and angular momentum vectors 
of k -th orbit in the constellation, respectively. 




FIG. 7 RIGID ROTATION MATRIX WITH FREE SELECTION OF 
AXIS AND ANGLE OF ROTATION 

Node vector and eccentricity vector of each orbit in the 
constellation can be expressed in terms of orbit 
elements as follows. 
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(18) 



where i k , Q. k , a> k are the inclination, the right ascension 

of the ascending node and the argument of perigee of 
k -th orbit, respectively. Then, considering rotating- 
plane lattice flower constellations with circular orbits 
and same orbit inclination, we obtain 



e k =f 2 (i,n k ) 

Rigid rotation matrix is given by 

R(a, tf>) = I 3k3 cos <j> + (1 - cos <j>) aa T + [ax] sin ( 



(19) 



(20) 



where a and <f> are axis and angle of rotation, 
respectively [Schaub et al. 2009]. 

Using the rigid rotation matrix, angular momentum 
and eccentricity vectors can be obtained as follows 



h; = R(M) h k = gl a' k ,n' k ) 

e' k =R(M) e k =g 2 (i k ,,n k .) 



(21) 



Results 

Optimal Radio Occultation Constellations 

1) String-of-pearls Polar Constellation 

The first constellation analyzed for radio 
occultation mission is a single-orbit, conventional, 
string-of-pearls constellation whose design 
parameters are given in Table 2. There is no single 
unique optimal solution because inclination i is 
independent of the active time when there is one 
orbital plane in constellation. In order to provide 
global coverage the orbit inclination is set to 90° 
(polar orbit). 

TABLE 2 LFC PARAMETERS OF STRING-OF-PEARLS POLAR 
CONSTELLATION 



Parameters 


Optimal Values 


[N ,N m ] 


[1,12] 


N c 




T 


1.51 (hr) 


i 


90 (deg) 



Fig. 8 illustrates the satellite trajectories of the 
string-of-pearls polar constellation as seen from the 
Earth-Centered Inertial (EG) frame. All satellites of 
this constellation have same orbital plane and all 
two consecutive satellites see each other, as shown 
in the figure. This polar orbit constellation provides 



global coverage due to the rotation of the Earth. 
However, long time is required to complete the 
global coverage and the coverage density is 
strongly latitudinal dependent with maximum 
observation at poles. 



ECI cecrdinsca 




FIG. 8 STRING-OF-PEARLS POLAR CONSTELLATION 

2) Plane Lattice Flower Constellations 

Although all satellites of a constellation don't have 
same orbital plane, it is possible to design 
constellations whose satellites lie on the same plane 
using different orbital planes. Compared with 
conventional string-of-pearls polar constellation, 
these constellations have advantage of relatively 
short time to cover global area although they can 
not provide global coverage. The disadvantage is 
that the deployment of this constellation may 
require more than a single launch. These 
constellations are called plane lattice flower 
constellations. The LFC parameters of designed 
plane lattice flower constellations are provided in 
Table 3. First one is called rotating-plane lattice flower 
constellation because axis of plane is coning about a 
fixed axis. The other is called fixed-plane lattice 
flower constellation because the plane where the 
satellites do not change its orientation. 

TABLE 3 LFC PARAMETERS OF PLANE LATTICE FLOWER 
CONSTELLATION 



Parameters 


Rotating-plane 
LFC 


Fixed-plane 
LFC 


[N o ,N sij ] 


[6,2] 


[12,1] 


K 


1 





T 


1.51 (hr) 


1.51 (hr) 


i 


21.78 (deg) 


76.94 (deg) 



Figs. 9-10 show the satellite trajectories of plane 
lattice flower constellations as seen from the ECI 
frame. 
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FIG. 9 ROTATING-PLANE LATTICE FLOWER CONSTELLATION 




FIG. 10 FIXED-PLANE LATTICE FLOWER CONSTELLATION 

When all satellites of a constellation lie in the same 
plane, then the rank of a matrix made of the 
position vectors must be 2. This position matrix is 
given by Eq. (22). 



M = [f, 



f 1 



(22) 



where wis the total number of satellites. Since the 
rank of a matrix is equal to the number of non-zero 
singular values, the singular values of a quadratic 
form of the matrix have been computed. 

Minimum singular values of the position matrix for 
various inclinations, [N o ,N m ,N c ] = [6,2,1] and 
r = 1.51(hr) are shown in Fig. 11 for the rotating- 
plane LFC. The rotating-plane LFC is really planar 
for low inclinations. Of course, for fixed-plane LFC, 
the rank of the position matrix is always 2, as we 
expected. 

Note that, using designed conditions, rotating- 
plane LFC can be used for low inclinations. In case 
of fixed-plane LFC, minimum singular values for 
various inclinations are always zero. 




FIG. 11 MINIMUM SINGULAR VALUES FOR VARIOUS 
INCLINATIONS, [Af,,iV M ] = [6,2], N c =l and T = 1.51 (hr) IN 
ROTATING-PLANE LFC 

Two Coverage Performance Evaluations 

1) First Coverage Estimation Result- sum of Areas 
of Simple and Complex Quadrilaterals 

T 

In this simulation, we set a time step to . 

r 100 

Single planar swath by two satellite of rotating- 
plane lattice flower constellation with single time 
step is displayed in Fig. 12. Quarter period 
coverage can be estimated by summing all planar 
swaths by all satellite combinations during a 
quarter period. However, the estimated coverage is 
not net coverage because there are overlapping 
areas as shown in Fig. 13. 



% 
.: 

% -30 





A > 




> 
























J *r 















Longitude (degree) 

FIG. 12 SINGLE PLANAR SWATH BY TWO SATELLITES WITH 
SINGLE TIME STEP 




Numbei of coverage 

FIG. 13 A QUARTER PERIOD COVERAGE OF ROTATING-PLANE 
LATTICE FLOWER CONSTELLATION 
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2) Second Coverage Estimation Result- Using 
Uniform Distributed Points on a Sphere 

In order to verify method using uniform 
distribution of point on a sphere, total coverage are 
estimated using various number of uniform 
distributed points on a sphere. Single period 
coverage percent and errors obtained with different 
number of points on string-of-pearls polar 
constellation are in Table 4. Because single planar 
swath of string-of-pearls polar constellation is 
narrower than other constellations, two coverage 
estimation methods have most significant 
difference in this constellation. Therefore, two 
coverage estimation methods are applied to string- 
of-pearls polar constellation. Errors are calculated 
by comparing with result which is obtained by sum 
of quadrilaterals (8.64%). 

TABLE 4 COMPARISON OF TOTAL COVERAGE OBTAINED 
WITH DIFFERENT NUMBER OF POINTS 



Number of Points 


Total Coverage (%) 


Error (%) 


2,560 


15.24 


76.4 


5,120 


11.56 


33.8 


10,240 


9.65 


11.7 


20,480 


8.96 


3.7 


81,920 


8.86 


2.6 


163,840 


8.67 


0.3 


1,310,720 


8.65 


0.1 



From these results, we can conclude that method 
using uniform distribution of point on a sphere 
needs appropriate number of points to give 
accurate results. Simulation with 1,310,720 points 
was performed in order to estimate net coverage. 
Single period simulation results of all designed 
constellations are shown in Figs. 14-16. Results by 
uniform distribution of point on a sphere are 
shown in Earth-Centered Earth-Fixed (ECEF) frame 
and have good agreement with 2D-view results of 
first coverage estimation method. 



3) Comparison of Coverage Performances 

As mentioned above, net coverage can be obtained 
by method with uniform distribution of points in 
sphere. Comparison of coverage performances by 
designed constellations are provided in Table 
5.Plane lattice flower constellations provides global 
coverage with relatively short time while string-of- 
pearls polar constellation require long time to cover 
global area. As shown in Fig. 15 and Fig. 16, 
rotating-plane lattice flower constellation provide 
good performance in a low latitude region while 
fixed-plane lattice flower constellation is suitable in 
a high latitude region. 

TABLE 5 COMPARISON OF COVERAGE PERFORMANCES (NET 
COVERAGE, %) 



Time 


String-of- 
pearls Polar 


Rotating-plane 
LFC 


Fixed-plane 
LFC 


1 period 


8.63 


36.20 


77.27 


0.5 day 


68.95 


36.90 


87.87 


1 day 


90.32 


36.92 


88.09 



Re-oriented Rotating-Plane Lattice Flower 
Constellations 

Using the relationship Eq. 21, we obtain orbit elements 
of the re-oriented constellation.Note that, re-oriented 
constellation has different inclination while original 
plane lattice flower constellation has identical 
inclination. 

Figs. 17 and 18 show the satellite trajectories and 
single period coverage of the re-oriented plane lattice 
flower constellations. Note that the re-oriented plane 
lattice flower constellations provide service in specific 
regions. The re-oriented fixed-plane flower 
constellation provides good performance in mid 
latitudes areas including highly populated regions 
such as the United States, Europe and Northeast Asia. 
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FIG. 14 SINGLE PERIOD COVERAGE OF STRING-OF-PEARLS POLAR CONSTELLATION 
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FIG. 15 SINGLE PERIOD COVERAGE OF ROTATING-PLANE LATTICE FLOWER CONSTELLATION 
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FIG. 16 SINGLE PERIOD COVERAGE OF FIXED-PLANE LATTICE FLOWER CONSTELLATION 




(a) (b) 



FIG. 17 RE-ORIENTED ROTATING-PLANE LATTICE FLOWER CONSTELLATION (a) THE SATELLITE TRAJECTORIES (b) SINGLE PERIOD 

COVERAGE IN 2D- VIEW 




y<kmt mm ' Longitude (degree) 

(a) (b) 

FIG. 18 RE-ORIENTED FIXED-PLANE LATTICE FLOWER CONSTELLATION (a) THE SATELLITE TRAJECTORIES (b) SINGLE PERIOD 

COVERAGE IN 2D- VIEW 



88 



Frontiers in Aerospace Engineering Volume 2 Issue 2, May 2013 



www.fae-journal.org 



Conclusions 

In this work, we analyzed and designed constellations 
for radio occultation mission of the Earth atmosphere. 
The Genetic Algorithms technique with 2-D Lattice 
Flower Constellation theory has been used to find the 
best constellation for the mission. Three different 
satellite constellations are found by maximizing active 
time with minimum distance constraint. Optimal 
solutions have not only conventional string-of-pearls 
polar constellation but also novel plane lattice flower 
constellation which have good coverage characteristics. 
In order to evaluate coverage performance, geometries 
of linear and planar swaths have been investigated 
and total coverage has been estimated. To count 
overlapping areas, uniform distribution of points on a 
sphere is introduced and net coverage of three 
different satellite constellations have been estimated 
by checking observation conditions of all uniformly 
distributed points on the sphere. It has been found 
that the proposed plane lattice flower constellations 
provide good performance for time to cover global 
area although they can't provide coverage for entire 
Earth. The rotating-plane lattice flower constellation 
and the fixed-plane lattice flower constellation have 
good coverage characteristics in low and high latitude 
region, respectively. In addition, the re-oriented plane 
lattice flower constellations provide good coverage 
performance in specific regions. 
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